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Galaxy structures are certainly fractal up to a certain crossover scale Ao. A clear 
determination of such a scale is still missing. Usually the conceptual and practical 
implications of this property are neglected and the structures are only discussed 
in terms of their global amplitude. Here we present a compact summary of these 
implications. First, we discuss the problem of the identification of the crossover 
scale Ao and the proper characterization of the scaling. We then consider the im- 
plications of these properties with respect to various physical phenomena and to 
the corresponding characteristic values, i.e. ro, erg, Q, etc. These implications cru- 
cially depend on the value of Ao, but they are still important for a relatively small 
value, say Ao ~ 50ft. - 1 Mpc. Finally we consider the main theoretical consequences 
of these results. 



1 Introduction 

Nowadays there is a general agreement about the £atvt-that galactic structures 
are fractal up to a distance scale of ~ 50/i -1 AfpcE3i23 and the increasing in- 
teres^atj jttjt^ tJ|^J'^^t^d | versus homogeneous distribution of galaxy in the last 



yearliJ'H'cBQ'EJ'ESLlcS'ESl has focused, mainly on the determination of the homo- 
geneity scale Ao[] The main point in this discussion is that galaxy structures 
are fractal no matter what is the crossover scale, and this fact has never been 
properly appreciated. Clearly, qualitatively different implications are related 
to different values of Ao . 

• Characterization of scaling properties. 

Given a distribution of points, the first main question concerns the possi- 
bility of defining a physically meaningful average density. In fractal-like 
systems such a quantity depends on the size of the sample, and it does 
not represent a reference value, as in the case of an homogeneous dis- 
tribution. Basically a system cannot be homogeneous below the scale 
of the maximum void present in a given sample. However the complete 
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where all these materials have 
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statistical characterization of highly irregular structures is the objective 



The major problem from the point of view of data analysis is to use 
statistical methods which are able to properly characterize scale invari- 
ant distributions, and hence which are also suitable to characterize an 
eventual crossover to homogeneity. Our main contributiono@'C3, in this 
respect, has been to clarify that the usual statistical methods (correla- 
tion function, power spectrum, etc.) are based on the assumption of 
homogeneity and hence are not appropriate to test it. Instead, we have 
introduced and developed various statistical tools which are able to test 
whether a distribution is homogeneous or fractal, and to correctly charac- 
terize the scale-invariant properties. Such a discussion is clearly relevant 
also for the interpretation of the properties of artificial simulations. The 
agreement about the methods to be used for the analysis of future sur- 
veys such as the Sloan Digital Sky Survey (SDSS) and the two degrees 
Fields (2dF) is clearly a fundamental issue. 

Then, if and only if the average density is found to be not sample-size de- 
pendent, one may study the statistical properties of the fluctuations with 
respect to the average density itself. In this second case one can study 
basically two different length scales. The first one is the homogeneity 
scale (Ao), which defines the scale beyond which the density fluctuations 
become to have a small amplitude with respect to the average density 
(dp < p). The second scale is related to the typical length scale of the 
structures of the density fluctuations, and, according to the terminology 
used in statistical mechanics E3, it is called correlation length r c . Such 
a scale has nothing to do with the so-called " correlation length" used in 
cosmology and corresponding to the scale £( r o) = E3, which is instead 
related to Ao if such a scale exists. 

• Implication of the fractal structure up to scale Ao . 

The fact that galactic structures are fractal, no matter what is the ho- 
mogeneity scale Ao , has deep implication on the interpretation of several 
phenomena such as the luminosity bias, the mismatch galaxy-cluster, the 
determination of the average density, the separation of linear and non- 
linear scales, etc., and on the theoretical concepts used to study such 
properties. We discuss in detail some of these points. We then review 
some of the main consequences of the power law behavior of the galaxy 
number density, by relating various observational quantities (e.g. ro, erg, 
ft, etc.) to the length scale Aq. 
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We also note that the properties of dark matter are inferred from the ones 
of visible matter, and hence they are closely related. If now one observes 
different statistical properties for galaxies and clusters, this necessarily 
implies a change of perspective on the properties of dark matter. 

• Determination of the homogeneity scale Xq. 

This is, clearly, a very important point which is at the basis of the un- 
derstanding of galaxy structures and more generally of the cosmological 
problem. We distinguish here two different approaches: direct tests and 
indirect tests. By direct tests, we mean the determination of the condi- 
tional average density in three dimensional surveys, while with indirect 
tests we refer to other possible analyses, such as the interpretation of 
angular surveys, the number counts as a function of magnitude or of dis- 
tance or, in general, the study of non-average quantities, i.e. when the 
fractal dimension is estimated without making an average over different 
observes (or volumes). While in the first case one is able to have a clear 
and unambiguous answer from the data, in the second one is only able 
to make some weaker claims about the compatibility of the data with a 
fractal or a homogeneous distribution. For example the papers of Wu 
ct aEj and Nusser & Laha-\E3 mainly concern with compatibility argu- 
ments, rather than with direct tests. However, also in this second case, 
it is possible to understand some important properties of the data, and 
to clarify the role and the limits of some underlying assumptions which 
are often used without a critical perspective. We do not enter here in the 
details of the discussion about real data (see e.g. Joyce et al. 1999, Wu 
ct al. 1999), however in the last section we consider separately the case 
(i) A w bQh^Mpc, (ii) A ~ ZOOh^Mpc and (iii) A « lOOOh^Mpc, 
briefly discussing the main theoretical consequences. 

2 Characterization of scaling properties 

In this section we describe in detail the correlation properties of a fractal 
distribution of points having eventually a crossover towards homogeneity (see 
Gabrielli & Sylos Labini, 2000 for a more exhaustive discussion of the subject), 
as the distribution of galaxies is thought to be (see Sylos Labini, Montuori & 
Pietronero, 1998 and Wu, Lahav & Rees, 1999 for two opposite views on the 
matter). 



Let 




(1) 
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be the number density of points in the system (the index i runs over all the 
points) and let us suppose to have an infinite system. If the presence of an 
object at the point r\ influences the probability of finding another object at 
r*2, these two points are correlated. Hence there is a correlation at the scale 
distance r if 

G(r) = (n(O)n(r)) ^ (n) 2 (2) 

where we average over all occupied points of the system chosen as origin and 
on the total solid angle, supposing statistical isotropy. On the other hand, 
there is no correlation if 

G(r) = (nf. (3) 
2.1 Homogeneity scale and correlation length 

The proper definition of A , the homogeneity scale, is the length scale beyond 
which the average density becomes to be well-defined, i.e. there is a crossover 
towards homogeneity with a flattening of G(r). The length-scale Ao is related 
to the typical dimension of the largest voids in the system. On the other hand, 
the correlation length r c separates correlated regimes of the fluctuations with 
respect to the average density from uncorrelated ones, and it can be defined 
only i£ a crossover towards homogeneity is shown by the system, i.e. if Ao 
cxist£3. In other words r c defines the organization in geometrical structures 
of the fluctuations with respect to the average density. Clearly r c > Ao: only 
if the average density can be defined one may study the correlation length of 
the fluctuations from it. In the case Ao is finite and then (n) > 0, in order 
to study the correlation properties of the fluctuations around the average and 
then the behaviour of r c , we can introduce the correlation function 

_ W0Wr»- W ' . (4) 

In the case of a fractal distribution, the average density (n) in the infinite 
system is zero, then G(r) = and Ao — oo and consequently £(r) is not 
defined. In this case the only well defined quantity characterizing the two 
point correlations is the function r(r)ed€l: 

r(r) = hm <"W0k (5) 
Rs^oo {n)R B 

where R s is the size of the a generic finite sample of the system, {■■■) R indicates 
the average over all the points of the sample as origins, hence (n)_R s is the 
average density of the sample. This function measures the average density 
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of points at a distance r from another occupied—point, and this is the reason 
why it is called the conditional average density cil. Obviously in the case of a 
distribution for which Ao is finite T(r) provides the same information of G(r), 
i.e. it characterizes the correlation properties for r < Ao and the crossover to 
homogeneity. 

A very important point is represented by the kind of information about the 
correlation properties of the infinite system which can be extracted from the 
analysis of a finite sample of it. In Pictronero (1987) it is demonstrated that 
even in the super-correlated case of a fractal the estimate of T(r) extracted 
from a finite sample, is not dependent on the sample size R s , providing a good 
approximation of that of the whole system. Clearly this is true a part from 
statistical fluctuations due to the fact that in a finite sample the average over 
the all possible origins is an average over a finite number of points, while in 
the global infinite system the average is over an infinite number of points. In 
fact, r(r) extracted from a sample can be written in the following way: 

i N i ff+Ar 

r « = ivg4^i n(?< + (6) 

where N is the number of points in the sample, n(fi + f") is the number of 
points in the volume element d 3 r' around the point ?i + ¥ and Ar is the 
thickness of the shell at distance r from the point at f^. 

Therefore, from an operative point of view, having a finite sample of points 
(e.g. galaxy catalogs), the first analysis to be done concerns the determination 
of r(r) of the sample itself. Such a measurement is necessary to distinguish 
between the two cases: (1) a crossover towards homogeneity in the sample 
shown by a flattening of T(r), and hence an estimate of Ao < R s and (n); 
(2) a continuation of the fractal behavior. Obviously only in the case (1), 
it is physically meaningful to study the correlation function f(r) (Eq.Q), and 
extract from it the length scale ro (£( r o) = 1), which is related to the intrinsic 
homogeneity scale Ao- The functional behavior of £(r) with distance gives 
instead information on the correlation length of the density fluctuations. 

2.2 The case of a fractal distribution (R s <C AoJ 

Hereafter we study the three-dimensional case, i.e. d = 3, and we suppose that 
the sample is a sphere of radius R s . Obviously, this choice is not a restriction. 

Let us analyze the case R s <C Ao < r c . This is the so called "fractal" 
case, and it is compatible with both the situation of Ao finite, but R s <C Ao (a 
sample-size which is smaller than the homogeneity scale), or the situation in 
which Aq — > oo, i.e. the case of a fractal distribution at any scale. 
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It is simple to shoviSiS that in this case (and in a spherical sample), 
Eq.|| becomes 

r(r) = (7) 

with B = N/Rj. Note that B is independent on the sample size: in fact, 
by changing R s , N in average scales as iif . This shows the aforementioned 
assertion that T(r) is practically independent on the sample-size. On the other 
hand, it is possible to show that B is related approximatively to the average 
distance between nearest neighbors points in the system Ej: 

where r e is the Euler's gamma-function. 

2.3 The "standard" correlation function for a fractal distribution 

As already mentioned, in the fractal case {R s <C Ao), the sample estimate of 
the homogeneity scale, through the value of r for which the sample-dependent 
correlation function £(r) (given by Eq.^) is equaLto 1, is meaningless: This 
estimate is the so-called "correlation length" ro cB in the standard approach 
of cosmology. As we discuss below, ro has nothing to share with the true 
correlation length r c . Let us see why ro is unphysical in the case R s <C Ao. 
£(r) is given operatively by 

« r) = ™>^_ 1 = £M (9) 

( n )R s \ n 'Rs 

The basic point in the present discussio is that the mean density of the 
sample, (n.)fl s , used in the normalization of £(r), is not an intrinsic quantity 
of the system, but it is a function of the finite size R s of the sample. 

In fact, from EqJ^—Lhe expression of the £(r) of the sample in the case of 
fractal distributions isEJ 

D ' — 



« r > = ¥(«:) - 1 - < 10 > 

From Eq.^ it follows that ro is a linear function of the sample size R s 
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and hence it is a spurious quantity without physical meaning but it is simply 
related to the sample's finite size. 

We note that the amplitude of T(r) (Eq.0) is related to the lower cut-off 
of the fractal £ by Eq.[| while the amplitude of £(r) is related to the upper 
cut-off (sample size R s ) of the distribution. This crucial difference has never 
been appreciated appropriately. 

Finally we stress that in the standard analysis of galaxy catalogs the fractal 
dimension is estimated by fitting £(r) with a power law, which instead, as 
one can see from Eq.|lO|, it is power law only for r <C ro (or £ 3> 1). For 
larger distances there is a clear deviation from the power law behavior due 
to the definition of £(r). Again this deviation is due to the finite size of 
the observational sample and does not correspond to any real change in the 
correlation properties. It is easy to see that, if one estimates the exponent at 
distances r < ro, one systematically obtains a higher value of the correlation 
exponent due to the break of £(r) in a log-log plot. To illustrate more clearly 
this we compute the log derivative of Eq.|lO] with respect to log(r), indicating 
D — 3 with 7 and its estimate with 7': 

_ d(log(j(r)) 2rgr-T 
7 dlog(r) 2r^r-7-l 7 ' 1 j 

where ro is defined by Eq|0]. The tangent to £(r) at r = tq has a slope 7' = 27. 
This explain why itJias-heen found in galaxy and cluster catalogs that t~ 2 
by the £(r) analysisl3oc3 instead of 7 ~ 1 found with the T(r) analysisc3. 



2.4 The case of a fractal distribution with a crossover to homogeneity (R s ^> 

Let us now analyze the case of a fractal with a crossover to homogeneity. In 
Coleman & Pietronero (1992) a very simple approximation has been used to 
describe such a situation which we discuss in more detail below. n 

By defining (n) R = N/V with V — 47ri?j?/3, it is simple to seda that 
the behavior of T(r) in our sample is fractal (i.e. T(r) is a power law) up to 
a certain distance Ao, and then it flattens, becoming homogeneous at scales 
R s > r > A : 

r(r) = ^-r D - 3 for I < r < A 

(13) 

r(r) ~ (n) Rs for A < r < i? s , 

where is the estimation of the average density in the sample of size R s . 

That is, (n) R does not depend on r if Ao -C r < apart small amplitude 
fluctuations. In Eqjl^ the detailed approach to homogeneity depends on the 
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specific properties of the fluctuations around the average density, i.e. it is 
determined by r c . Hence, the statistical properties of the density fluctuations 
determine how good is the estimation of the average density throught (n) R . 
From the definition of the function £(r) we can finccJ 

sw *(£T'fe)- (I4) 

Note that the amplitude of £(r) is determined by the homogeneity scale Ao 
which has been previously extracted from T(r), and that in this approxima- 
tion ro ~ Ao . The function £(r) characterizes the correlations among the 
fluctuations of the distribution with respect to the average density. It is im- 
portant to clarify that these fluctuations must be both positive and negative, 
in fact the integral over the whole sample of Eq.|| must be 0. Hence the func- 
tion / must be oscillating, and in the case A < r C! R Sl it should present 

an exponential cut-off at r rs r c . We have that, when £(Ao) — 1, the density 
fluctuations begin to become small with respect to the average density (n), 
but if Ao < r < r c they are still well correlated among them. Only for r 3> r c 
the fluctuations are not correlated. 

Let us now consider the case Ao <C R s < t c . This situation is compatible 
with the following two situations: r c finite, but larger than i? s , and the case 
r c — > oo. In both cases, £(r) of our sample should be a power law modu- 
lated by an oscillating function g{r) which describes the positive and negative 
fluctuations with respect to the average density, 

^ r )={£) 1 g{r) - (15) 

In such a situation the (positive and negative) fluctuations from the average 
density are of all sizes and they do not have any intrinsic characteristic scale: 
this is a critical system (see Gaite et al., 1999 for a more detailed discussion). 
The only intrinsic scale of the system is then Ao, the length-scale beyond which 
r(r) flattens and the fluctuations are small with respect to the average. 

Let us suppose to be in the case in which T(r) flattens at a certain Ao <C R s - 
We can then evaluate the correlation function £(r) of the sample via Eq^. At 
this point we can clarify how to interpret the eventual cut-off shown by £(r). 



If the cut-off scale is well below R s , we can be sure that it is a good 
estimate of the intrinsic correlation length r c ; 



• if the cut-off is at a scale r ~ R s we can have two cases: it represents 
an "intrinsic" cut-off with r c ~ R s , or it is only a finite-size effect due 
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r (Mpc/h) 

Figure 1: The conditional average density T(r) for a distribution which has a power law 
behavior at small scales (r < Ao) , followed by a transition to homogeneity at the scale 
Ao = lO/i" 1 Mpc. The behavior of the flattening depends on the correlation properties of the 
density fluctuations, i.e. on the functional behavior of £(r). The dotted line corresponds to a 
system which has a finite correlation length r c = 30h~ Mpc, while the solid line describes a 
system whose density fluctuations present correlation over all scales. From the r(r)-analysis 
alone it is possible to compute Ao but not r c . 

to the fact that from Eq|| |£(i? s )| =0. In order to distinguish between 
these two possibilities, it is necessary to increase the sample size and to 
look at the behavior of the cut-off scale. If it increases proportionally 
to R Sl then it is a finite size effect. Otherwise if it does not change, it 
represents the estimate of the intrinsic correlation length r c . 

In Fig.|l| we show two possible behaviours of the flattening of T(r), while in 
Fig.^it is shown the corresponding £(r) (we neglect for simplicity the oscillating 
term which must be present, and we have considered the situation R s — > oo). 
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r (Mpc/h) 



Figure 2: The £(r) correlation function for the distributions shown in the previous figure. 
With this analysis it is possible to compute the correlation length r c (finite or infinite) of 
the density fluctuations. The dotted line corresponds to an exponential decay, and hence to 
a finite value of the correlation length (r c = 30/i _1 Mpc) while the solid line corresponds to 
an infinite correlation length, and hence to a power law behavior. The length scale at which 
£(»*o) = 1 gives a reliable estimation of Aq. 



10 



2.5 About the amplitude of £(r) 



We note that if A <C R s , \q has nothing to share with questions like "which 
is the typical size of structures in the system?" or "up to which length-scale 
the system is clusterised?" . The answer to this question is strictly related to 
r c and not to Ao- The length scale r c characterizes the distance over which 
two different points are correlated (clusterised p3. In fact, this property is not 
related to how large are the fluctuations with respect to the average (Ao), but 
to the length extension of their correlations (r c ). 

To be more specific, let us consider a fixed set of density fluctuations. 
They can be superimposed to different values of a uniform density background. 
The larger is this background the lower Ao, but obviously the length scale of 
the correlations (r c ) among these fluctuations is not changed, i.e. they are 



clusterised inder 
One can sec 



ently of the background. 

that a linear amplification of £(r) such that 

£'(r) = M(r) 



(16) 



doesn't change r c (which can be finite or infinite) but only Ao, i.e. if A > 1 we 
need larger subsamples to have a good estimation of (n) , but the characteristic 
length (correlation length) of the structures is not changed. 



2. 6 Homogeneity scale and the size of voids 

Basically Ao is related to the maximum size of voids: the average density 
will be constant, at least, on scales larger than the maximum void in a given 
sample. Several authors have approached this problem by looking at voids 
distribution. For example El- Ad and Piran (1997) have shown that the SSRS2 
and IRAS 1.2 Jy. redshift surveys are dominated by voids: they cover the 
~ 50% of the volume. Moreover the two samples show very similar properties 
even if the IRAS voids are ~ 33% larger than SSRS2 ones because they are not 
bounded by narrow angular limits as the SSRS2 voids. The voids have a scale 
of at least ~ 40 4- 50h~ 1 Mpc and the largest void in the SSRS2 sample has a 
diameter of ~ 60h~ 1 Mpc, i.e. comparable to the Bootes void. The problem 
is to understand whether such a scale has been fixed by the samples' volume, 
or whether there is a tendency not to find larger voids: in this case one would 
have a (weaker evidence) for the homogeneity scale. In any case, we note that 
the homogeneity scale cannot be smaller than the scale of the largest void 
found in these samples and that one has to be very careful when comparing 
the size of the voids to the effective depth of catalogs. For example in the 
Las Campanas Redshift Survey, even if it is possible to extract sub-samples 
limited at ~ 500h~ 1 Mpc, the volume of space investigated is not so large, as 
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the survey is made by thin slices. In such a situation a definitive answer to 
the dimension of the of voids, and hence to the existence of the homogeneity 
scale, is rather difficult and uncertain. 

2. 7 Luminosity Bias 

We would like to stress again that, even if the fractal behavior breaks at a 
certain scale Ao, the use of £(r) is in anyhow inconsistent at scales smaller 
than Ao- We illustrate below an example of the confusion due to the use of 
£(r) when r]UXo. 

From the use of the £(r) analysis, it has been found that r is different in 
different volume (hereafter VL) samples. In particular it has been foundcl that 
deeper is the VL sample, larger is the value of r$ . As the deeper VL samples 
contain brighter galaxies, this fact has been interpreted as a real physical 
phenomenon, leading to the idea that more brighter galaxies are more strongly 
clustered than fainter ones, in view of their larger correlation amplitude: this 
is the so-called luminosity segregation phenomenon PH 37 !. In other words, the 
fact that the giant galaxies are "more clustered" than the dwarf ones, i.e. 
that they are located in the peaks of the density field, has given rise to the 
proposition that larger objects may correlate up to larger length scales and 
that the amplitude of £(r) is larger for giants than for dwarfs one. The deeper 
VL samples contain galaxies which are in average brighter than those in the 
VL samples with smaller depths. As the brighter galaxies should have a larger 
correlation length the shift of ro in different samples can be related, at least 
partially, with the phenomenon of luminosity segregation. 

As previously discussed there are two problems with such a model: (i) 
The amplitude of £(r) in an homogeneous distribution, does not give any in- 
formation about the clustering " strength" . It is instead related to the local 
amplitude of the fluctuations with respect to the average density, (ii) The 
amplitude of £(r) has a physical meaning only in the case Ao is found to be 
finite and smaller than the sample's size. This is clearly not the case up, at 
least, to ~ 50h- 1 Mpc. 

A natural explanation of the scaling of ro is then the fractal behavior 
of galaxy distribution, and more specifically the fact that ro is a fraction of 
the sample's size in the fractal case. The fact that giant elliptical galaxies 
are located in the core of rich clusters, and other morphological properties 
of this kind, .can be naturally related to the multifractal properties of matter 
distributioil§c3. In such a case, bright galaxies are more strongly clustered 
than fainter ones in view of the fact that their fractal dimension is smalle£3. 
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2.8 Power Spectrum of density fluctuations 



The problems with the standard correlation analysis also show that the prop- 
erties of fractal correlations have not been really appreciated. These problems 
are actually" far more serious and fundamental than mentioned, for example, 
by LandyEj and the idea that they can be solved by simply taking the Fourier 
transform is once ipqpa a proof of the superficiality of the discussion. We have 
extensively shown oc3 that the power spectrum of the density fluctuations 
has the same kind of problems which £(r) has, because it is normalized to 
the average density as well. The density contrast 8(r) = 8p{r)/(p) is not a 
physical quantity unless the average density is demonstrated to exist. More 
specifically, like in the case of £(?"), the power spectrum (Fourier Transform of 
the correlation function) is affected by finite size effects at large scale: even 
for a fractal distribution the power spectrum has not a power law behavior 
but it shows a large scale (small k) cut-off which is due to the finiteness of the 
sampleE-l Hence the eventual detection of the turnover of the power spectrum, 
which is expected in CDM-like models to match the galaxy clustering to the 
anisotropics of the CMBR, must be considered a finite size effect, unless a clear 
determination of the average density in the same sample has been done. 

Essentially all the currently elaborated models of galaxy formation El as- 
sume large scale homogeneity and predict that the galaxy power spectrum 
(hereafter PS), which is the PS of the density contrast, decreases both toward 
small scales and toward large scales, with a turnaround somewhere in the 
middle, at a scale A/ that can be taken as separating "small" from "large" 
scales. Because of the homogeneity assumption, the PS amplitude should be 
independent on the survey scale, any residual variation being attributed to 
luminosity bias (or to the fact that the survey scale has not yet reached the 
homogeneity scale). However, the crucial clue to this picture, the firm de- 
termination of the scale A/, is still missing, although some surveys do indeed 
produce a turnamund scale around 100 /i -1 AfpctrE2l. Recently, the CfA2 sur- 
vey analyzed byE3 (hereafter PVGH) (and confirmed by SSRS2liil - hereafter 
DVGHP), showed a n = —2 slope up to ~ 30/i _1 Mpc, a milder n ~ — 1 slope 
up to 200 h~ l Mpc, and some tentative indication of flattening on even larger 
scales. PVGH also find that deeper subsamples have higher power amplitude, 
i.e. that the amplitude scales with the sample depth. 

In the following we argue that both features, bending and scaling, are a 
manifestation of the finiteness of the survey volume, and that they cannot be 
interpreted as the convergence to homogeneity, nor to a PS flattening. The 
systematic effect of the survey finite size is in fact to suppress power at large 
scale, mimicking a real flattening. Clearly, this effect occurs whenever galaxies 
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have not a correlation scale much larger than thepSiixvey size, and it has often 
been studied in the context of standard scenarioscJo. We push this argument 
further, by showing that even a fractal distribution of matter, which never 
reaches homogeneity, shows a sharp flattening and then a turnaround. Such 
features are naxtially corrected, but not quite eliminated, when the correction 
proposed by t£l is applied to the data. We show also how the amplitude of 
the PS depends on the survey size as long as the system shows long-range 
correlations. 

The standard PS (SPS) measures directly the contributions of different 
scales to the galaxy density contrast Sp/p. It is clear that the density contrast, 
and all the quantities based on it, is meaningful only when one can define a 
constant density, i.e. reliably identify the sample density with the average 
density of all the Universe. In other words in the SPS analysis one assumes 
that the survey volume is large enough to contain a homogeneous sample. When 
this is not true, and we argue that is indeed an incorrect assumption in all the 
cases investigated so far, a false interpretation of the results may occur, since 
both the shape and the amplitude of the PS (or correlation function) depend 
on the survey size. 

Let us recall the basic notation of the PS analysis. Following Peebles 
we imagine that the Universe is periodic in a volume V u , with V u much larger 
than the (presumed) maximum homogeneity scale. The survey volume V G V u 
contains N galaxies at positions rj, and the galaxy density contrast is 

Sir) = ^ - 1 (17) 
n 

where it is assumed that exists a well defined constant density n, obtained 
averaging over a sufficiently large scale. The density function n{r) can be 
described by a sum of delta functions: n(r) = Y2iLi S^(r — fl) . Expanding 
the density contrast in its Fourier components we have 

jeV 

where 

W(k) = 1 / drW{r)e li? (19) 

is the Fourier transform of the survey window W{r), defined to be unity inside 
the survey region, and zero outside. If £(r) is the correlation function of the 
galaxies (£(f) =< n(r)n(0) > /h 2 — 1), the true PS P(k) is defined as the 
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Fourier conjugate of the correlation function £(r). Because of isotropy the PS 
can be simplified to 

P{k) = 4tt / Z( r f-^j^r 2 dr . (20) 



The variance of 5? is 

<\h\ 2 >=l + ^m- (2i) 

The first term is the usual additional shot noise term while the second is the 
true PS convoluted with a window function which describe the geometry of 
the sample (PVGH) 



p (k) = I < M 2 > \W(k-k')\ 2 d 3 k'. (22) 



with 



V 

P{k) = I dk'P{k')F{k - k') , (23) 

F(k-k') = ^\W{k-k')\ 2 . (24) 

We apply now this standard analysis to a fractal distribution. We recall the 
expression of £(r) in this case is 

C(r) = [(3- 7 )/3](r/i2 a )-T-l, (25) 

where 7 = 3 — D. A key point of our discussion is that that on scales larger 
that R s the £(r) cannot be calculated without making assumptions on the 
distribution outside the sampling volume. 

As we have already mentioned, in a fractal quantities like £(r) are scale 
dependent: in particular both the amplitude and the shape of £(r) depend the 
survey size. It is clear that the same kind of finite size effects are also present 
when computing the SPS, so that it is very dangerous to identify real physical 
features induced from the SPS analysis without first a firm determination of 
the homogeneity scale. 

The SPS for a fractal distribution model described by Eqj2^ inside a sphere 
of radius i?, is 



, f,fls sin(fer) 
P(k)= / 4tt- 



kr 



3- 7 fj_y _ { 

R, 



2 , a k (R s )R 3 - D b k (R s ) 
rdr = fc5 fci-. 

(26) 
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Notice that the integral has to be evaluated inside R s because we want to 
compare P{k) with its estimate in a finite size spherical survey of scale R s . 
In the general case, we must deconvolve the window contribution from P{k); 
R s is then a characteristic window scale. Eq.p6| shows the two scale-dependent 
features of the PS. First, the amplitude of the PS depends on the sample depth. 
Secondly, the shape of the PS is characterized by two scaling regimes: the first 
one, at high wavenumbers, is related to the fractal dimension of the distribution 
in real space, while the second one arises only because of the finiteness of the 
sample. In the case of D = 2 in Eq.|26| one has: 

47T 

a k (R s ) = —(2 + cos(kR s )) (27) 

and 

b k (R s ) =47rsm(fciZ s ). (28) 

The PS is then a power-law with exponent —2 at high wavenumbers, it flattens 
at low wavenumbers and reaches a maximum at k w 4.3/i? s , i.e. at a scale A rj 
1.45i? s . The scale at which the transition occurs is thus related to the sample 
depth. In a real survey, things are complicated by the window function, so that 
the flattening (and the turnaround) scale can only be determined numerically. 

In practice one has several complications. First, the survey in general is not 
spherical. This introduces a coupling with the survey window which is not easy 
to model analytically. For instance, we found that windows of small angular 
opening shift to smaller scales the PS turnaround. This is analogous to what 
happens with the correlation function of a fractal: when it is calculated in small 
angle surveys, the correlation length r decreases. Second, the observations are 
in redshift space, rather than in real space. The peculiar velocities generally 
make steeper the PS slope E3 with respect to the real space. Third, in a fractal 
the intrinsically high level of fluctuations makes hard a precise comparison 
with the theory when the fractal under study is composed of a relatively small 
number of points. 

3 Implications for cosmology 

We now consider some implications for cosmology of the scaling properties of 
galaxy distribution, up to a lenght scale Xq. For example, we consider more 
specifically Ao ~ 50h~ 1 Mpc. 

3. 1 Estimation of the average luminosity and mass density 

From the studies of Large Scale Structures (LSS) of galaxies and galaxy clus- 
ters, one would like to estimate the average density of visible matter and then 



1G 



to infer the one of the whole (visible plus dark, i.e. all the matter in chistcrised 
objects) matter distribution. While for the first we have direct estimations, for 
the second we have only indirect methods, especially at large scales, based on 
some assumptions, which can be tested by looking at the distribution of what 
is observable, i.e. visible matter. 

We briefly describe how to do such a measurement in galaxy redshift cata- 
logs directly, from the knowledge of galaxy positions and luminosity (for a more 
detailed discussion see Sylos Labini, 2000). In this case, and by measuring the 
Mass-to-Luminosity ratio, one can infer, from the average luminosity density, 
the average mass density. The new point we address more specifically is that 
galaxies are fractally distributed up to a certain crpasoyer scale Ao. As there is 
still some controversy about the value of Ao oe3u'HE-3 we give the estimation 
as a function of Ao- We stress that the way this estimation is performed, is 
substantially different from the usual onec3, because in such a case the frac- 
tal behavior is not considered at all, and one assumes a perfect homogeneous 
distribution at relatively small scale (Ao ~ 5 -j- Wh Mpc). This situation 
is clearly not the one corresponding to the more "optimistic" estimation of 
the homogeneity scale Ao E3e3'E9. The other assumption usually made-is that 
galaxy positions are independent on their luminosity. We have showicl that, 
although such an assumption cannot describe local morphological properties 
of galaxy distribution^, it works rather well in the available galaxy redshift 
surveys. 

The estimation of the the average density we are able to make depends 
hence on two parameters. The first one is the homogeneity scale Ao and the 
second is the Mass-to-Luminosity ratio. We can give.-an upper limit to £1 by 
taking the highest {M. / C) c — 300/i observed up to nowEJ (in clusters of galaxies) 
to be universal across all the scales, and by considering a lower limit for the 
homogeneity scale Ao = 50h~ 1 Mpc. We compute the critical (M/C) cr it, i-e. 
the Mass-to-Luminosity ratio needed to have fi = 1. As the others, also this 
parameter depends on the homogeneity scale Ao- 

3.2 Average luminosity density from galaxy catalogs 
Let 

(v(r,L))dLd 3 r = (t>{L)(T{r))dLd z r = Ar D - 3 L a e~"ft d 3 rdL (29) 

be the average number of galaxies in the volume element d 3 r at distance r from 
a observer located on a galaxy, and with luminosity in the range [L, L + dL\. 
In Eq.^9| we have used the fact that the galaxy luminosity function has been 
observed to have the so-called Schechter shape with parameters L* (luminos- 
ity cut-off) and a (power law index) which can be determined experimentally. 
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The conditional average space density (r(r)) has a power law behavior corre- 
sponding to a fractal dimension D (which eventually can be a function of scale, 
and hence can approach to D = 3 at a scale Ao)- Both the fractal dimension 
D and the overall amplitude A can be determined in redshift surveys. Hence 
(f(r, L)) is a function of four parameters: L*, a, D, A. Moreover we note that 
by writing (i/(r, L)) as a product of the space density and of the luminosity 
function we have implicitly assumed that galaxy positions are independent on 
galaxy luminosity. 

We would like to estimate the average luminosity density in a sphere of 
radius R and volume V(R) placed around a galaxy, and defined as 

i pR poo / t> \ D—3 

which is R dependent as long as the space density shows power law behavior 
(i.e. LL< 3). By considering M* = 19.53 (i.e. L* = 1.0 • 10 10 /r 2 Z©), a = 
— 1.05E3 and by estimating the prefactor A (Ea.p9[) and the fractal dimension 
(D sa 2) in galaxy redshift samples we obtaiio 

^ (10) w 2 ■ 10 8 hL Q /Mpc 3 . (31) 

We now estimate the density parameter in terms of the critical density, 

Pc = 2.78 • 10 n /i 2 M Q /Mpc 3 (32) 

where M Q is the solar mass. By considering the product of the mass-to- 
luminosity ratio (in solar and h units) and the average luminosity density 
given by Eqs. [3C)|]3l"| , we obtain 

W = (6±2).10-^ (^) , (33) 

where Ao is the scale where the crossover to homogeneity occurs (it can also be 
Ao = oo, and in such a case O(oo) = 0). Note that in view of the dependence 
of A4/C on h, Eq.|33| does not depend on the Hubble's constant. 

Let us pow suppose that M./C ~ lOh as it has been derived by Faber and 
GallengherEy. We obtain 

n(Ao)«6.10-»( i A_)" 1 . (34) 

If galaxy distribution turns out to be homogeneous at AopSi 10h~ 1 Mpc then 
fl w 6 ■ 10^ 3 as it is obtained in the standard treatment c3. If, instead, the 
crossover to homogeneity lies at 100h~ 1 Mpc, we obtain w 6 ■ 10 -4 . 
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From Eq.[33| and if galaxy distribution turns out to be homogeneous at Ao , 
we obtain that the critical Mass-to-Luminosity ratio (such that fi(Ao) = 1) is 
given by 

so that if Ao = Wh 1 Mpc one obtains. {M./£) c rit ~ 1600ft. (which is again con- 
sistent with the usual adopted valuecS) while if Ao = 100h~ 1 Mpc {M./ '£) C rit ~ 
16000ft, which is about two orders of magnitude larger than the highest A4/C 
observed in clusterised objects. For an intermediate value of Ao = h0h~ 1 Mpc 
one obtains (M/C) cr it ~ 8000ft. 

For what concerns the analysis of galaxy clusters it is often used a value 



(M/C) c ~ 300/iBEl which, by using EqM gives 



fi(Ao)«2.10- 1 (^ I ^_) . (36) 

Such an estimation is based on the fact that the Mass-to-Luminosity ratio 
found in clusters is representative of all the field galaxies. This is a very strong 
assumption: this means that (M/C) g = (M/C) c which is not supported by 
any observation. The usually adopted value of Q = 0.2u can be derived from 
Eq|36| by assuming Ao = 10h~ 1 Alpc. In the case the crossover to homogeneity 
occurs at Ao = 100h~ 1 Mpc we have that fi(Ao) = 0.02 if we consider (M/C) c ~ 
300/i to be "universal". 

Under the assumptions: 

(i) galaxies are homogeneously distributed at scales larger than Ao « 
5h~ 1 Mpc, 

(ii) the A4/C of galaxies is the same of clusters, that is galaxies should 
contain a factor ~ 30 more-dark matter than what it is observed with the 
study of the rotation curvescil, 

(Hi) Galaxy positions are independent of galaxv-Juminosity: such a as- 
sumption is not strictly valid, but it has been testecc3 to hold rather well in 
the available samples, 

we get the following upper limit to fL If we assume that A4/C — 300h 
across all the scales, and that Ao ~ 50/i _1 Afpc, which we consider to be a lower 
limit for the homogeneity scale, we get from Eq.[36| 

n(50Mpc/h, 300hM Q /L Q ) < 0.04 . (37) 



The direct estimation with galaxies (Eqs. p3[p4| ) gives lower a value, the 
reason being the strong assumption of taking A4/C = 300ft. as representative 
of field galaxies (see FigJ|). 
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Figure 3: We show the behavior of f2(r, A4/C) and Q c (t) (direct estimation from the Mass 
and number density of galaxy clusters), in the case h = 1, for three different values of 
M/C = 1, 12,300. If A = b0h~ 1 Mpc we get an upper limit Q, = 0.04. In the figure the 
value of the Hubble's constant is set h = 1. (From Sylos Labini, 2000) 
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3.3 Homogeneity scale and primordial nucleosynthesis constraints 

Let us now discuss the previous results in relation with the nucleosynthesis 
constrain Q BBN = 0.015ft, -2 , deriving in such a w.aY an upper limit for Ao 
which is consistent with such a scenario. We may se£3 that if Ao ~ 100h 2 Mpc 
then f2(Ao, (-M/£)c) ~ flf BN ■ Such a fact has two important implications. 

First of all one does not need the presence of non baryonic dark matter to 
reconcile local observations of matter contained in galaxies and galaxy clusters 
(fliocai) with the primordial nucleosynthesis constraint. This is an important 
point as the existence of non baryonic dark-matter has been inferred also, but 
not only, by the discrepancy Qi oca i/Q BBN a, which, we show, is not observed 
if Ao > l00h 2 Mpc. There is still place, in this picture, for a uniform back- 
ground of non-baryonic dark matter, which has completely different clustering 
properties from the ones of visible matter. 

The second point concerns the baryon fraction in clusters. According to the 
usual root, in order to-.be consistent with U BBN , one has to force = 0.1 -=-0.3 
from clusters analysis!!. Here we can argue as follows. By assuming a very 
large M/C ~ 300ft, and that it is universal across all scales, and by assuming 
that all dark matter is made of baryons, if the crossover to homogeneity occurs 
at scales larger than ~ 150h 2 Mpc then there is not enough baryonic matter 
to satisfy the nucleosynthesis constraint. The situation gets clearly worst if 
one takes into account that not all the mass in clusters is baryonic or that 
M/C < 300ft, i.e. it is not considered to be universal across all scales, or 
that Ao 3> 100ft 2 Mpc. There are then two different possibilities: the first is to 
study the case of a non homogeneous primordial nucleosynthesis which could 
lower the limit on il BBN , given the observed abundances of light elements, 
and the second would be to have a more uniform background of baryonic 
dark matter. This seems rather unlikely because it would have very different 
clustering properties with respect to visible mass and it is very difficult to find 
a dynamical explanations for such a segregation. 

3.4 Where does linear approximation hold ? 

In the standard picture, the properties of dark matter on cosmological relevant 
scales r > 5h~ 1 Mpc, are inferred from the observed properties of visible matter 
and radiation. Now one studies change in these properties, i.e. the presence 
of fractal co rre lations, and in this respect they will have consequences on dark 
matter toaJjtj. For example, the determination of the mass density including 
dark matter has been performed on the basis of the linear theorylij. Here the 
problem is: beyond which scale can linear theory be considered as a useful 
approximation ? 
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In other words, the dynamical estimates use gravitational effects, of depar- 
ture from a strictly homogeneous distribution on the motion of objects such as 
galaxies considered as a test particle. A completely different situation occurs 
if Ao is larger than the scales at which linear approximation is..usually adopted. 
For example, methods based on the Cosmic Virial TheoremEJ, distortions in 
redsthft surveys Ej, local group dynamics c2l, the "reconstructed" peculiar ve- 
locity field from the density field (e.g POTENT-like methodO), are clearly not 
useful at scales r <C Ao. Up to now it is implicitly assumed Ao ~ 5-i-10h Mpc 
and all these methods are considered to be valid on larger scales. If, for exam- 
ple Ao ~ 50h~ 1 Mpc then it is not possible to interpret peculiar velocities in-thc 
range ls-30h~ 1 Mpc through the linear approximation (as it is usually don(e3), 
unless there is a background of dark matter which is uniform beyond a certain 
small scale. In .such a situation however the estimates would be differenta from 
the usual onesc2l. 

Let us review some simple relations between the conditional average den- 
sity r(r) and the usual correlation function £(r). If T(r) has a power law 
behavior up to a scale Ao and then it presents a crossover to homogeneity, it is 
simple to show that (in the case of D — 2) the length scale at which £(ro) = 1 
is of the order of 

ro « £ (38) 

where the_ exact relation between ro and Ao depends on the details of the 
crossover E3. However Eq.|3^ gives a reasonable order of magnitude in the 
general case. In such a situation it is also simple to compute <j(r, Ao) = (N(r) — 
(N)/(N)), i.e the_amplitude of the fluctuation with respect to the average at 
the scale r <C AoC3 

*(r,Ao)«^. (39) 

It has been found in various nearby surveys that a(8h~ 1 Mpc) = 1. However 
if the crossover to homogeneity occurs at Ao we have that 

aiSh-HlpcXo)-^ (40) 
16 

For example if Ao ~ 15h~ 1 Mpc then ro w 5/i _1 Mpc and a(8h~ 1 Mpc, 15h~ 1 Mpc) = 
1; otherwise if A « 50h" 1 Mpc then r « IQh^Mpc and a(8h^ 1 Mpc, Ibh^Mpc) 
3. 

Linear approximation holds in the linear regime, when the amplitude of 
the density contrast is small, i.e. for a <C 1. We have that a — 1 at a scale 
rv=i(Ao) 

r CT=1 (A ) « ^ (41) 
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For instance, if the crossover to homogeneity occurs at Ao = hQh~ 1 Mpc one has 
that cr(Ao) = 1 at r a= i — 25h~ 1 Mpc. In such a situation all the estimationsjQf 
the density parameter at smaller scales based on the linear approximation^, 
and hence based on the untested assumption of linearity, are not correct. 

We have studied in detaicJ the gravitational force distribution in a fractal 
structure. Its behaviour can be understood as the sum of two parts, a local 
or 'nearest neighbours' piece due to the smallest cluster (characterised by the 
lower cut-off A in the fractal) and a component coming from the mass in other 
clusters. The latter is bounded above by the scalar sum of the forces 

(\F\) < lim [ L GPm } r) Avr 2 dr~L D - 2 (42) 
L ^°° J a r 2 

so that for D < 2 it is convergent, while for D > 2 it may diverge. If there is 
a divergence, it is due to the presence of angular fluctuations at large scales, 
described by the three-point correlation properties of the fractal. For the dif- 
ference in the force between two points the local contribution will be irrelevant 
well beyond the scale A, while it is easy to see that the 'far-away' contribution 
will now converge as L D ~ 3 , and its being non-zero is a result of the absence 
of perfect sphcricaLsymmetry. We have then applied such a result to the case 
of an open univers£H in order to compute the expected deviations from a pure 
linear Hubble flow. 

4 Discussion and Conclusions 

We now present a short discussion about the perspective of our work: the 
cosmological implications of the fractal behavior of visible matter crucially 
depend on the crossover scale Ao, but, no matter what is the actual value of 
such a scale, we have some important consequences from the theoretical point 
of view. We may identify three different scenarios. 

• (1) The fractal extends only up to ~ 30 -f- \QQh~ 1 Mpc. Thia is the min- 
imal concept which begins to be absorbed in the literatures but, some- 
times, without considering its real consequences. The standard approach 
to galaxy distribution has identified very small " correlation lengths" , 
namely 5h~ 1 Mpc for galaxies and 25h~ 1 Mpc for clustery These num- 
bers (which were supposed to know with high precisioro) are anyhow 
inconsistent with the fractal extending a factor 2-^4 more. We have 
shown that this inconsistency is conceptual and not due to incomplete 
data or week statistics. Hence, in this hypothesis one has to abandon all 
the concepts related to these length scales. These are: 
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(i) The estimate of the matter density in clusterised objects (visible + 
dark), which has been claimed to be fl s» 0.2 -=- 0.3, decreases by one 
order of magnitude or more. 

(ii) The normalization of N-body simulations is usually performed to 
some length-scale or amplitude of fluctuations, which are related to 
hbT x Mpc and 25h~ 1 Mpc. 

(iii) Concepts like the galaxy-cluster mismatch and the related luminosity 
bias, as well as the understanding of the clustering via the bias parameter 
b (i.e. linear or non-linear - "stochastic bias" - amplification of £(r)) loose 
any physical meaning. 

(iv) The interpretation of the velocity field is also based on the linear 
approximation which cannot certainly hold at scales smaller than 30 4- 
50h- 1 Mpc. 

(v) The reconstruction of the three dimensional properties from the an- 
gular data suffers of the same untested assumption of homogeneity. 

In summary major modifications are necessary for the origin and dynam- 
ics of large scale structures and for the role of dark matter. However the 
structures may still formed via gravitational instability, in the sense that 
they are not necessarily primeval. 

• (2) The fractal extends up to 300-^500/i -1 Afpc. In this case the standard 
picture of gravitationally induced structures after the electromagnetic 
decoupling is untenable. There is no time to create such large scale 
correlated structures via gravitational instability, starting from Gaussian 
initial conditions. More string consequences are clearly important for 
what concerns the amount of matter in clusterised objects. 

• (3) The fractal extends up to > 1000h~ 1 Mpc, and homogeneity does not 
exist, at least for what concerns galaxies. In this extreme case a new 
picture for the global metricfii is then necessary. 

For some questions the fractal structure leads to a radically new perspec- 
tive and this is hard to accept. But it is based on the best data and analyses 
available. It is neither a conjecture nor a model, it is a fact. The theoretical 
problem is that there is no dynamical theory to explain how such a fractal Uni- 
verse could have arisen from the pretty smooth initial state we know existed 
in the big bang. However this is a different question. The fact that something 
can be hard to explain theoretically has nothing to do with whether it is true 
or not. Facing a hard problem is far more interesting than hiding it under 
the rug by an inconsistent procedure. For example some interesting attempts 
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to understand why gravitational clustering generaLes-scale-invariant structures 
have been recently proposed by de Vega et aLSESE-l Indeed this will be the 
key point to understand in the future, but first we should agree on how these 
new 3d data should be analyzed. In addition, the eventual crossover to homo- 
geneity has also to be found with our approach. If for example homogeneity 
would really be found say at ~ lQ0h~ 1 Mpc, then clearly all our criticism to 
the previous methods and results still holds fully. In summary the standard 
method cannot be used neither to disprove homogeneity, nor to prove it. One 
has simply to change methods. 
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